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Abstract 

The enclosure method was originally introduced for inverse problems of concern- 
ing non destructive evaluation governed by elliptic equations. It was developed as 
one of useful approach in inverse problems and applied for various equations. In 
this article, an application of the enclosure method to an inverse initial boundary 
value problem for a parabolic equation with a discontinuous coefficients is given. A 
simple method to extract the depth of unknown inclusions in a heat conductive body 
from a single set of the temperature and heat flux on the boundary observed over 
a finite time interval is introduced. Other related results with infinitely many data 
are also reported. One of them gives the minimum radius of the open ball centered 
at a given point that contains the inclusions. The formula for the minimum radius 
is newly discovered. 
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1 Introduction 



Assume that we have a set of the pair of the temperature field on the boundary of a 
heat conductive body and the corresponding heat flux across the boundary of the body 
over a finite time interval. A part in the body that has a different conductivity from the 
known reference one is called an inclusion. In this paper we consider the problem: what 
information about inclusions in the body can one extract from the set? The solution 
to this problem may have possible application to non destructive evaluation by thermal 
imaging. We study this problem from a mathematical point of view and aim at seeking 
an analytical approach for extracting information about the location and shape of the 
inclusions. 
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Let Q be a bounded domain of with a smooth boundary. We denote the unit 
outward normal vectors to by the symbol v. Let T be an arbitrary fixed positive 
number. 

Given / = f{x,t), {x,t) G dQx]0, T[ let u = u{x,t) be the solution of the initial 
boundary value problem for the parabolic equation: 

Ut-W ■ -fVu = inQx ]0, T[, 
7VM-t/ = /onai]x]0, T[, (1.1) 
u{x, 0) = in Jl, 

where 7 = 7(0;) = (7ij(a;)) satisfies 

(Gl) for each i,j = 1, 2, 3 -fij{x) = -fji{x) e L°°(Q); 

(G2) there exists a positive constant C such that 7(0;)^ • ^ > C|^p for all ^ G and a. 
e. X &VL. 

This paper is concerned with the extraction of information about "discontinuity" of 7 
from u and 7VM ■ v on (9nx]0, T[ for some / and an arbitrary fixed T < 00. However, we 
do not consider completely general 7. Instead we assume that there exists an open set D 
with a smooth boundary such that D C Vt and 7(0;) a.e. x E Vl\D coincides with the 
3x3 identity matrix Is and satisfies one of the following two conditions: 

(Al) there exists a positive constant C such that —(7(2^) — ^3)^ • C ^ for all ^ G R^ 

and a.e. x & D; 

(A2) there exists a positive constant C such that {j{x) — I^)^ ■ ^ > C"|^p for all ^ G R'^ 
and a.e. x & D. 

Write h{x) — j{x) — 1^ a.e. x G D. In this paper we consider the following problem: 

Inverse Problem. Assume that both D and h are unknown. Extract information about 
the location and shape of D from a set of the pair of temperature ^(a;, t) and heat flux 
f{x,t) for {x,t) G x]0, T[. 

The D is a model of the union of unknown inclusions where the heat conductivity 
is anisotropic, different from that of the surrounding homogeneous isotropic conductive 
medium. The problem is a mathematical formulation of a typical inverse problem in 
thermal imaging. 

Elayyan-lsakov [4] investigated the uniqueness issue of this type of problem. As a 
corollary of their uniqueness theorem we know that the lateral Neumann-to Dirichlet map: 
/ I — )■ M|9nx]o,r[ uniquely determines D together with h inside D ii VL \ D is connected 
and h is given by 6/3 with a smooth function h on D. However, their purpose is to recover 
the full information about the location and shape of D and for the purpose their proof 
requires infinitely many pairs of the temperature and heat flux on SJlxjO, T[ even just 
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for determining a single point on dD. This shows a difficulty of obtaining the detailed 
image of inclusions from boundary measurements. Note also that in [1, 2] an approach 
to Inverse Problem in a onc-spacc dimensional case which is based on the idea of the 
prohe method introduced by Ikehata [6] has been proposed. However, the procedure is 
quite complicated compared with the original probe method and needs huge number of 
measurements. 

In this paper we mainly seek a simpler method that yields a partial or rough information 
about the location and shape of D from u{x, t) on dVlx ]0, T[ for a single fixed heat flux or 
explicit heat fluxes prescribed on 917 x ]0, T[. We think that this type of study provides 
us with a knowledge about good heat fluxes on the boundary of the body to get such 
an information. In [11] we have already developed an argument based on the enclosure 
method which was originally introduced for elliptic equations in [7, 8] to derive two types 
of formulae in the case when the inclusion has the zero conductivity^ that is a cavity. The 
argument yields the values of the support function of the cavity at a given direction and 
the distance of a given point outside the body to the cavity from the temperature fields 
and special explicit heat fluxes. In this paper, we will see that the argument also works 
for the inclusion case and yields also a new information: the minimum radius of the open 
ball centered at a given point that contains the inclusions. 

The main new point of this paper is: an introduction of another argument which is 
also based on the enclosure method and gives a formula which has not been considered in 
[11]. It makes use of a single set of a general heat flux and the corresponding temperature 
field on the surface of the body over a finite time interval. It yields a depth of unknown 
inclusions in a heat conductive body from the surface of the body. Wc do not prescribe 
any explicit heat flux on the surface of the body, instead assume a standing behaviour for 
a given heat flux. 

Note that in Theorem 2.1 of [9] the enclosure method has been applied to a one-space 
dimensional version of Inverse Problem. Therein complex exponential solutions of the 
backward heat equation with a large parameter are used. In this paper we use only real 
exponential solutions. 

1.1 A formula with a general heat flux 

The new point of this paper is a derivation of the following formula which can be considered 
as the main result of this paper. It makes use of a single set of a heat flux and the 
corresponding temperature on dVLx ]0, T[ and gives a pre-knowledge about the location 
of inclusions. 

Theorem 1.1. Let f = f{x,t) satisfy: there exists a /i. G R such that 

rT rT 

< inf liminf r'^ / e-^*/(a;, ^)c^^ < sup limsupr'' / e'^^ f(x,t)dt < oo (1.2) 
xedci T — >-oo Jo xedn t — *-oo Jq 

and that the function 

fT 

dfl3x\ — >g{x;T)= / e~''*f(x,t)dt 

Jo 

is continuous. 
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Let u — Uf{x,t) be the weak solution of (1-1) for this f and let v = Vg{x;T) be the 
solution of 

(A-r)v = infl, 

av 

— = q(x: t) ondfl. 

Then, there exists a tq > such that 

• if (Al) is satisfied, then for all r > Tq 



I I 

J da Jo 



e""^* {vg{x; T)f{x, t) — Uf{x, t)g{x; r)) dtdS < 0; 
• if (A 2) is satisfied, then for all r > tq 

/ / e~''^ {vg{x;T)f{x,t) - Uf{x,t)g{x;T))dtdS > 0. 

Jon Jo 



I an Jo 

In both cases the formula 

rT 



Am ^ log 



r — >oo 



an Jo 



e 



-Tt 



Vg{x; T)f{x, t) — Uf{x, t)g{x; r)) dtdS 



is valid, where 

dist{D, dQ) = mi{\y - x\\y e dQ, x e D}. 



= -dist{D,dn), 

(1.4) 



Note that Varadhan [14] considered the asymptotic behaviour as r — > oo of the 
solution of the problem 

(A-T)v = OinQ, 
V — 1 ondfl. 

He used the behaviour to establish the short time asymtotics of the heat kernel. See 
also [13] and references therein for the subject itself. Theorem 1.1 shows that this type of 
solutions can be applied to inverse initial boundary value problems for parabolic equations 
over a finite time interval. 

In [10] Ikehata considered an inverse obstacle scattering problem whose governing 
equation is given by the wave equation in three dimensions. The observation data are 
given by a wave field measured on a known surface surrounding unknown obstacles over a 
finite time interval. The wave is generated by an initial data with compact support outside 
the surface. Applying the idea of the enclosure method, he estabhshed an extraction 
formula of the distance from a given point outside the surface to obstacles from the data. 
To establish the formula he made use of the solution v e H^(R^) of the inhomogeneous 
modified Helmholtz equation 

(A-T^)v + f(x) = inR^ 

where f{x) is an initial data of the wave field. Thus the equations in (1.3) correspond to 
this equation. However, in contrast to the solution of this equation, that of (1.3) has not 
an explicit form in general. In this paper, we solve (1.3) by using the potential theory 
and study its behaviour as r — > oo to get a necessary estimate. 
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1.2 Other three formulae with special heat fluxes 



If one uses special heat fluxes, then one can explicitly obtain more information about the 
location and shape of D. The idea for the derivation of the following formulae come from 
[11]. 

The second result is the following. 

Theorem 1.2. Given u e S'^ let f be the function of {x,t) e 9Qx]0, T[ having a 
parameter r > defined by the equation 



. , , dv . . , . 
f{x,t;T) = —{x;r)cp{t), 



;i.5) 



where v{x; r) = e"^^'^ and e -^'^(0, T) satisfying the condition: there exists e R such 
that 



lim inf 



/ e-^^cp(t)dt 
Jo 



> 0. 



(1.6) 



Let Uf = Uf{x,t) be the weak solution of (1.1) for f — f{x, t; r) and hoiuj) — sup^.^^, x-u;. 
Then the formula 



lim — ;= 



is valid. 



log 



Jg^J^ ^ (^{x;T)f{x,t;T) -Uf{x,t)^{x;T)^ dtdS 



hoicj), (1.7) 



Note that: if if{t) is smooth on [0, T'[ with < T' < T and t = is not a zero point with 
infinite order of (fi{t), then (1.6) is satisfied for an appropriate fj, > 0. 

Next we choose a third solution of the equation (A — t)^; = in Q: given p e \ Q 

g-x/r|a;-p| 

v(x;t) = — — , X e Q. 

\x — p\ 

Using this v, we obtain the third formula. 

Theorem 1.3. Let p G R^ \ fl and replace v of f in (1.5) with the v above. Let 
Uf = Uf{x,t) be the weak solution of (1.1) for this f = f{x,t;T,p). Then assuming (1.6), 
one has the formula 



lim — -= 

r— >oo 2^T 



log 



an Jo 



-Tt 



' dv \ 

v{x; T)f{x, t; T,p) — Uf{x, t) — {x; r) J dtdS 



-dnip), 



where doip) denotes the distance from p to D, 

doip) =inf{|y-p| \ y &D}. 



Finally we introduce another formula which is also new and not given in [11]. Let 
1/ G R'^ be an arbitrary fixed point. We choose the function v given by 



and v{x;t) = 2\Jt at x = y. Note that the v{x;t) is smooth as the function of x and 
satisfies the modified Helmholtz equation in the whole space. Hence we can choose the 
reference point y E without any restriction. Note that Theorem 1.3 gives doip)', 
however, we have to take p e \ f2. 

Theorem 1.4. Let y G R"' and replace v of f in (1-5) with the v above. Let u 
he the weak solution of (1.1) for f 
formula 



f = Uf{x,t) 

f{x,t;T,y). Then assuming (1.6), one has the 



hm — j= 

T — >oo 2^T 



log 



an Jo 



-Tt 



v{x; T)f(x, t; T, y) - Uf(x, t)^(x; r) j dtdS 



Rniy), (1.9) 



where Rniy) — sup^g^ \x — y\. 

The theorem above makes use of a smooth solution of the modified Helmholtz equation 
that grows every points as r — > oo. The function Ru{y), y E VL is a. new comer and 
gives the minimum radius of the ball centered at y that contains D. Moreover we have 
the estimate of D from above as 



D c r\y^-R?{x eB?\\x-y\< Roiy)}- 



1.3 Construction of the paper 

A brief outhne of this paper is as follows. Theorem 1.1 is proved in Subsection 2.3 after 

formulating the notion of the weak solution of (1.1) together with a related estimate in 
subsection 2.1. The proof is based on an integral identity which is described in subsection 
2.2. Using the identity, we give an asymptotic representation formula of the integral 

lonlo'''^' (^<^'^^)f(^^t) -Uf{x,t)^{x;T)^ dtdS 

whose leading term is given by using two Neumann-to-Dirichlet maps for the operators 
A — r and V-7V — r in Q. Then with a help of a system of integral inequalities [10] which 
is widely used in previous applications of the enclosure method to elliptic equations [8] we 
see that the problem is reduced to giving some asymptotic estimates for the integral of 
the gradient of Vg over D. In some sense, this is an mdzrect verification of the hypothesis: 
Vg{x; r) ~ e^^/T'd-oni^) j- — y qq_ 'jj-^q estimates are stated in subsection 2.3 and their 
proof is given in subsection 2.4. It is based on the integral representation of Vg with a 
single layer potential over dQ. The proof of Theorems 1.2, 1.3 and 1.4 can be done along 
with the same line as [11] in which the case when dD is perfectly insulated is considered. 
For reader's convenience we describe an outline of the proof in section 3. In Appendix we 
give detailed proofs of four claims used in subsection 2.4. 



2 Extracting depth 

2.1 Preliminaries about the direct problem 

In this subsection, following [3] we describe what we mean by the solution (1.1). The 
presentation here is almost parallel to subsection 2.1 in [11]. 
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We put W{Q, T; H\n), {H^Q))') = {u e L'^{0, T; H\n)) \ u e L^O, T; (if^(Q))')}. 
Given / e L'^{0, T; H-^/^{dn)) we say that u e W{0, T; H^{n), {H^Q))') satisfy 

dtU-V ■ 'jVu = inf^xjO, T[, 

(2.1) 

-fWu-u ^ f ondnx]0, T[ 

in the weak sense if the u satisfies 

< u'{t), ip> + [ ^{x)Vu{x, t) ■ V<p{x)dx =< f{t), <p\9n > in (0, T), (2.2) 
Jn 

in the sense of distribution on (0, T) for all Lp G H^{Q) and a.c. t g]0, T[. We see that 
every u G W{0, T; H^{fl), {H^^fl))') is almost everywhere equal to a continuous function 
of [0, T] in L^{n) (Theorem 1 on p.473 in [3]). Further, we have: 

W{0, T;H\n), (H'in))') ^ C°([0, T];L\n)), (2.3) 

the space C°([0, T]; L^(f2)) being equipped with the norm of uniform convergence. Thus 
one can consider u{0) and u{T) as elements of L^(0). Then we see that given uq G L^(Q) 
there exists a unique m such that u satisfies (2.1) in the weak sense and satisfies the initial 
condition u{0) = Uq (Theorems 1 and 2 on p. 512 in [3]). 

Let Uq = 0. Remark 2 on p. 512 and Theorem 3 on p. 520 in [3] yields the continuity of 
u on. f: there exists a Ct > independent of / such that 

||ii||L2(o,T;i?i(n)) < C'T||/||L2(o,r;i?-i/2(an))- (2-4) 
Moreover, from (2.2) and (2.4) we have 

||'"'||L2(0,r;i?l(n)') < C'T||/||L2(o,T;i?-l/2(aQ))- 

This together with (2.3) and (2.4) yields one of the important estimates in the enclosure 
method: 

||'«(r)IU2(f2) < C'T||/||L2(o,r;i?-i/2(an)). (2.5) 

In the following subsection we denote by the weak solution of (2.1) with u{0) — and 
this is the meaning of the weak solution of (1.1). 



2.2 A basic identity 

Define 

rT 

Wf{x;T) — / e '^^Uf{x,t)dt, x & fl 

J 

and ^ 

gAx]T)= I e~''*f{x,t)dt, x E dQ, 
Jo 

where r > is a parameter. This type of transform has been used in the study [9] for the 
corresponding problem in a one-space dimensional case. 

In this subsection we derive an identity that connects the data for the parabolic equa- 
tion with the Cauchy data of the solutions of the modified Helmholtz type equations. 
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The Wf — w satisfies 

(V-7V-t)w = e-^^Uf{x,T) infl, 

^Vw ■ v — gf on 9Q. 

Let V = v{x) satisfy (A — r)^ = in ^2. Integration by parts yields 

/ lgfV-Wf^)dS=f{^-l3)Vv-Wwfdx + e~^^[uf{x,T)v{x)dx. (2.6) 

Let pf — p he the unique solution of the boundary value problem: 

(V • 7V - t)p = in Q, 

7Vp • v — gf on dfl. 

Set ef — Wf — Pf. Since we have 

L (^^" " = - • 

from (2.6) it follows that 



(2.7) 



+ / (7 — Is) ■ Ve/dx + e '^'^ / Uf{x,T)v{x)dx. 
Note that e/ = e satisfies 

(V • 7V - T)e = e~^^Uf{x, T) in Q, 
7Ve • = on 90. 



(2.8) 



Let Ris{t) and R^{t) denote the Neumann-to-Dirichlet maps on dO, for the operators 
A — r and V • 7V — r in respectively. We have 



(dv \ 

Rhi^) (^l^^^j =^lao, R^{T)gf=pf\ 



an- 



Since both Risir) and R^{t) are symmetric, we obtain from (2.7) 

L (^^" " ^^1^) ^ L ^^^^^^^^^ " ^^^^^^ (1^'^") 

+ / {j - Iz)Vv ■Vefdx + e''"'^ j Uf{x,T)v{x)dx. 



(2.9) 



This is our basic identity. In the proof of Theorems 1. 1 to 1.4 we show that, in some sense, 
one can ignore the second and third terms of this right-hand side. Thus the basic identity 
provides us a relationship between the boundary data for the parabolic equation over a 
finite time interval and the Cauchy data for the modified Helmholtz type equations. 
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2.3 Proof of Theorem 1.1 

Since the e satisfies (2.8), one gets 

L2(f2), (2-10) 

where C is a positive constant. Since / is independent of r, it follows from (2.5) and 
(2.10) that ||Ve||L2(fi) = 0(6"^^) as r — > oo. 

Now substitute v — Vg into (2.9). Prom (1.2) and (1.3) one gets ||t'g( • ; ''')||i/i(n) — 
0{t^^) as r — )■ oo. Prom these one gets the estimate on the second and third term in 
(2.9) as r — > oo: 

jja- h)^Vg ■ Vtfdx + e-^^ "/(^' T)vgdx = ©(r-'^e"^'^). 

Note that this is a very rough estimate, however, for our purpose it is enough; at this step 
we never make use of the assumption that 7(2:) = Is outside D. 
Summing up, we have obtained the asymptotic formula: 

/ {gVg-Wfg)dS^ I g{Rj,{r)-R,{r))gdS + 0{r-^e-^''). (2.11) 

JdCl JdCl 

The following system of inequalities is quite useful to give an estimation of the first term 
of this right-hand side. 

Proposition 2.1. Let ^0 and ^ satisfy (Gl) and (G2). Let u solve 

V • ^Vu — Tu — Q inD,, 



and V 



7 ■ v — g on dQ 

V • 7oVi' — TV — infl, 
jqVv ■ p = g ondfl. 



Then it holds that 



I {lo^ -T^)lQ^v --ioVvdx < f g{v-u)dS < [ {-f--fo)Vv-Vvdx. (2.12) 

Por the proof see [5]. In the present situation ^o{x) = I3 and ^{x) — I3 a.e. x e Q\D 
and thus from (2.12) we obtain 

/ {h-7-')^Vg-Vvgdx< [ g{Ri,{T) - R^{T))gdS < [ - h)Vvg -Vvgdx. (2.13) 

JD JdCl Jd 

Here we describe a key lemma whose proof is given in the next subsection. 
Lemma 2.1. There exist real numbers Ai and A2 independent of t such that 

limsupr^^e^^'^^^^^^'^") / \Wvg\''dx < 00 (2.14) 

T ^-OO J D 
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and 



liminfr^^e'^^^*^^^'^") / iVvJ^dx > 0. (2.15) 

T ^OO Jd 



From the proof one can choose Xi = 2fi — 1 and A2 = 2yU + 5/2, however, the exact 
values of Ai, A2 are not important for the derivation of formula (1.4) itself. 
Prom (2.11), (2.13) and (2.14) one gets 

limsupr^ie^^^i^^(^'^") / {gvg - Wfg) dS < 00. (2.16) 

Now consider the case when (Al) is satisfied. It follows from the right half of (2.13) 
and (2.15) that 

lirninfr^^e'^^^^^^-^'^") J^^g{Ri,{r) - R^{T))gdS^ > 0. 
This together with (2.11) gives 

liminf r^2g2V?dist(D,an) ^_ j^^ ^^^^ _ ^^^^ ^ q ^2.17) 

This imphes also that there exists a tq > such that for all r > tq 

- (gVg - Wfg) dS > 0. 

Next consider the case when (A2) is satisfied. Since — ^(x)^^ = ^(x)^^^'^{^j(x) — 
/3)7(a;)"^/^, one can find a positive constant C such that, for all ^ G R'^ (-^3~7(2^)^^)C'C ^ 
C|^|^. Hence a similar argument yields that 

lirninf r^2^^V^dist(D,dn) ^^^^ _ ^^^^ ^ q (2.18) 

and this implies that there exists a Tq > such that for all r > Tq 

j^^{gVg-Wfg) dS > 0. 
Now formula (1.4) is a direct consequence of (2.16), (2.17), (2.18) and the identity 

□ 

2.4 Proof of Lemma 2.1 

Let 11 be the constant in (1.2). Set v{x;t) = T^Vg{x;T) and 

g(x;T) = r e-^'f{x,t)dt, X e dn. 
Jo 

It suffices to prove (2.14) and (2.15) for v instead of Vg. 
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The V satisfies (A — t)v — m Q and dv/dv — g on dQ. In what follows we simply 
write V and g as v and g, respectively. We think that this makes no confusion. Thus from 
(1.2) one has: there exists positive constants C and Tq independent of x e such that, 
for all X e dfl and all r > Tq 

C-^ <g{x;T)<C. (2.19) 
Using the potential theory (cf. [12]), one has the expression 



1 r e~^l^~^l 

vix;T) = — -i){y]T)dSy, xeVt, 

27r Jdvi \x — y\ 



where -jt) G C{dfl) is the unique solution of the integral equation of the second kind 
on dQ: 

It is well known that the operator 

C{dn) 3ip^ S9n{T)ip e C{dn), 

where 

1 f d /e"^l^"^'l\ , 

Sdn{r)(p{y) = — ^ jr- ^{y')dSy., y e dQ 

2tt Jan ovy \ \y — y \ j 

is bounded and its operator norm has a bound 0(r~^/^) as r — )■ oo. Thus it follows 
from (2.19) and (2.20) that ■?/'(■ ;r) also has: there exists positive constants C and tq 
independent of a; G dVt such that, for all x G dVt and all r > Tq 

C"-^ < V'(x;t) < C". (2.21) 

Since 

e-V^\^-y\ ( 1 \ x-y 

\v\;x]T) = -— IWVy.r; 
zTT Jon 

one has 



Vv{x;t) = -— ij{y;T)- rK^+i 1 i ^"^^v^ 

ZTT Jdn \x — y\ \ \x~y\J \x~y\ 



[ \Wv(x;r)\^dx= [ dsJ rfS,. / rfa;e-^l^-^le-^l^-«'l$(x,y,y';T), (2-22) 
Jd Jan Jdn Jd 



where 

^(x,y,y';r) 



(2.23) 



(27r)2 \x — y\\x — y'\ \ \x — y\) \ \x — y'\J \x — y\ \x — y' 



From (2.21), (2.22) and (2.23) we can easily obtain (2.14) with Ai = -1. 

The problem is the proof of (2.15). We divide the integrand of (2.22) into two parts. 
Set do — dist (D, dQ) and M — {{x, y) E D x dn\\x — y\ — do}- It is easy to see that M 
coincides with the set of all {x, y) G dD x such that \x — y\ = do- 

In what follows we denote by Br{z) the open ball centered at a point z with radius 
R. Given S > define 

W5 = U^^o,yo)eM(D n Bs{xo)) x {OQ n Bs{yo)) x {OQ n Bs{yo)). 
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The set Ws is open in D x dD, x dD, and contains the set of all {x, y, y) with (x, y) ^ M.. 

Here we state two claims concerning the Ws whose proof is given in Appendix. 

Claim 1. Given e > there exists a 5i > such that for all {x, y, y') e Ws-^ it holds that 

X — y X — y' /I , 

1 r ■ 1 77 ^ 1 ~ ^5 \x — y\ < do + e, \x — y \ < do + e. 

\x — y\ \x — y'\ 

Claim 2. Given 6i > there exists a ^2 > such that if (x, y, y') e D x dQ x dQ\ Ws^, 
then \x — y \ + \x — y'\ > 2do + 62- 

Thus giving e = 1/2 in claim 1 and choosing the corresponding 61, we have: if 
{x,y,y') e Wsi, then 

X — y X — y' ^1 
\x — y\ \x — y'\ ~ 2 

It follows from this together with the left half of (2.21) and (2.23) that there exist constants 
Ci and To > such that, for all {x, y, y') e and r > tq 

$(x,y,y';T) > Cit. 

On the other hand, using the right half of (2.21), it is easy to see that: there exist positive 
constants C2 and ti > tq such that, for all {x, y, y') e D x dQ x dQ and r > Ti 

^x,y,y';T) < C2T. 

Now choose 82 in claim 2 corresponding to 5i already chosen. 
Then we have e-^d^-^l+l^"^'') < e'^^'^"'^^^ for any {x,y,y') e D x dft x dfl \Ws^. 
Hence dividing the integral (2.22) into Wsi and its compliment, one gets as r — > 00 

[ \Vv{x;T)\''dx>Ci[ dSydSy>dxe-^^^''-y^+^''-y'^^ + 0{e-^^'^-^^'). (2.24) 

Jd JWsi 

Choose a {xq, yo) E M.. It follows from the definition of Ws^ and the inequality \x — y\ + 
\x - y'\ < 2\x - xo\ + 2do + \yo - y\ + \yo - y'\ that 

/ dSydSy>dXe-^^^''-y^ + ^''-y'^^ 

> f_ dx f dSy f rfS'j./e-^d^-^/l+la.-s/'l) 

JddBs^ (xo) JdnnBs^ (yo) JdClDBs^ (yo) 

> g-2v/7do f e-2^l^-^°l(ix I / e-^^y-y^dSy 

JDnSg^ixo) yJdnnBsj^iyo) 

Now (2.15) with A2 = 5/2 is a direct consequence of this together with (2.24) and the 
following two claims. 
Claim 3. For all 5 > we have 

liminf(v^)=^ / e-^^l^-^oldx > 0. 

Claim 4. For all 5 > we have 

liminf r / e'^^y^-^USy > 0. 

r— >oo JdanBeiyo) 

For the proof of these claims see Appendix. 
□ 
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3 Outline of the proof of Theorems 1.2, 1.3 and 1.4 

Since from (1.5) we have 

Ov /"'^ 

it follows from (2.9) that 

^ ^ ^ ^ (3.1 

+ / {l - h)^v ■Vefdx + e"''^ I Uf{x,T)v{x;T)dx. 



By virtue of (2.5) one knows that for both v in Theorems 1.2 and 1.3 there exists a 
constant k such that \\uf{ ■ , ;T)\\l2^q) = 0{e'^^) as r — )■ oo. This together (2.10) yields 
that, as T — > oo 

/ (7 - /3)V^; • Vsfdx + e-^^ I Uf{x, T)v{x; T)dx = 0(6"^^/^). (3.2) 
Here we recall the following lemma. 

Lemma 3.1. There exist real numbers /^i, /X2, /^s G R sttc/i ^/ia^, /or all uj G 5*^, p G R^\f2 
and 1/ G 

liminfr'^ie-^^'^^^'^) / e^^^""dx > 0, (3.3) 

r— J-oo 7^, 

liminfr'^^e^^'^^^) / e-^^l^-^'lrfx > (3.4) 

and 

liminfr'^^e-^^-^^^^) / e^^l^-^lrfx > 0. (3.5) 

r — )-oo 

We do not mind the precise values of fii, fi2, fJ'?,- Every case can be reduced to the 
case when D is given by an open ball since we are assuming that dD is smooth. See 
[8] for the proof of (3.3) and [11] (or [10]) for the proof (3.4) and (3.5). Now it is a due 
course to see that a combination of (1.6), (2.12), (3.1), (3.2) and (3.3)/(3.4)/(3.5) yields 
(1.7)/(1.8)/(1.9). 
□ 



4 Conclusion and open problems 

We showed how the enclosure method can be applied to inverse initial boundary value 
problems over a finite time interval for parabolic equations with discontinuous coefficients. 
We established four types of formulae. It should be emphasized that in the all formulae 
the initial temperature field inside the body is assumed to be a known constant. We think 
that this is a natural condition and can be realized without special care in practice. In 
fact, just make it cold by using a refrigerator if the size of the body is not so large! 

Two of them are new in idea and yield: (I) a depth of unknown inclusions in a heat 
conductive body from the surface of the body with a single set of a heat fiux and the 
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corresponding temperature on the surface over a finite time interval (see Tfieorem 1.1); 
(II) the minimum radius of the open ball centered at a given point that contains unknown 
inclusions with a special explicit heat flux with a large parameter (sec Theorem 1.4). 

The point is the choice of the heat flux / and a solution v of the modified Helmholtz 
equation (A — t)v = in f2 in the integral 

In (I) first / is given and we choose v by solving the Neumann problem for the modified 
Helmholtz equation in Q whose Neumann data can be calculated from /; in (II) using a 
special v which is growing everywhere as r — > oo, we specify the form of /. 

The procedure suggested from (I) of extracting dist {D, dfl) is extremely simple and 
summarized as follows. 

(i) Give the heat flux / satisfying (1.2) for a e R across dQ over the time interval ]0, T[. 

(ii) Measure the temperature Uf{x,t) on 9fi over the time interval ]0, T[. 

(iii) Fix a large r > and compute the solution Vg of (1.3). 

(iv) Compute the quantity 

/ / e~'^^ \vix]T)fix^t) — utix^t)^{x-,T)\ dtdS 
Jon Jo \ Of J 

as an approximation of dist (D, dVL). 

If D is near surface and isolated in a small part, the information dist {D,dVt) may 
not be so useful, however, if D is deep inside or occupies a large part, then the set of all 
X eVL such that dist (D, dO) < ddQ{x) may give a good estimation of D from above. 

The method can be applied also to more complicated situations, for example, inclusions 
in a body with a known inhomogeneous isotropic or anisotropic conductivity apart from 
some technical difficulties or similar problems with acoustic /elastic/ electromagnetic waves, 
etc.. Such applications belong to our future study. 

A next challenging problem is: to clarify what information about D can be extracted 
from the asymptotic behaviour of integral (4.1) as r — > oo if / is fixed] the v is one of 
the three special solutions of the modified Helmholtz equation in Theorems 1.2, 1.3 and 
1.4, that is one of 

v{x]T,uj) = e^^"", X e R^ a; e 5^ 

q-s/t\x-p\ _ 

v{x;t,p) = -, X eR^\ {p}, p e \ n. 



gV^lx-i;! _ g-v^|a:-)/| 

v{x; r, y) = —— , x e R^ \ {y}, v{y; r, y) = 2^/^, y e Rl 

p y\ 

This remains open at the present time. 
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5 Appendix. Proof of claims 

5.1 Claim 1 

Define 

F(x,y,y')^(l-^^^-^^]+{\x-y\-do)+(\^-y'\-do), (x,y,y') eDxdQxdn. 
\ \x-y\ \x-y'\J 

Since |a; — > do and (x — y)/\x — y\ ■ {x — y')/\x — y'\ <1 for (x, y, y') E D x dfl x dfl, 
it suffices to prove that: given e > there exists a 5i > such that F{x,y,y') < e for 
all {x,y,y') G Wsi- Assume that this is not true. There exists a eo > and a sequence 
{xi,yi,yi) e Wi/i, I = 1,2, •• • such that F{xi,yi,y[) > eo. By the definition of Wi/i we 
know that, for each I there exists a {pi, qi) e A4 such that \xi — pi\ < l/l, \yi — qi\ < l/l 
and \y'i — qi\ < l/l. Since D and dfl are compact, one can choose a subsequence li, h- --- 

of / = 1, 2, ■ ■ ■ in such a way that the limits limj ^.oo xi- = x E D, limj ^ooPi = p E D, 

Mmj ^ooVij = y ^ dQ, limj ^ooVi^ = u' ^ dQ and lim^ ^oo % = Q ^ dQ exist. Clearly 

it holds that x = p and y = y' = q. Since At is closed, one gets {x, y) E A4 and thus 
\x — y\ — do. This together with y — y' gives F{x,y,y') — 0. On the other hand, since 
F{xi.,yi.,y'i.) — one gets F{x,y,y') > eo. This is a contradiction. 
□ 

5.2 Claim 2 

Assume that the statement is not true. There exist a Sq > and a sequence {xi,yi,y'i) € 
D X ^Vtx^V^\Ws^■,, / = 1, 2, ■ ■ ■ such that \xi~yi\ + \xi~y[\ < 2do + l/l. Since D and 90 are 
compact, if necessary replacing the sequence with a suitable subsequence, one may assume 

that the limits lim; xi = x E D, lim; 5.00 yi = y E dO, and lim^ 5.00 y'l = u' ^ dQ exist. 

Since Ws^ is open, one has (x, y, y') E D x dQ x dQ\ Ws^. 

On the other hand, since {\xi — yi\ — do) + {\xi — y'i\ — do) < l/l, \xi — yi\ > do and 
\xi — y'll > do, we obtain \x — y\ = \x — y'\ = do- This means that {x,y) E M. and 
{x, y') E Ai. Since x E dD, using local coordinates at x and y, one can easily show that 
{y — x)/do = and similarly {y' — x)/do = Ux- This yields y = y' and thus {x, y, y') E W^q. 
This is a contradiction. 
□ 

5.3 Claim 3 

Given (5 > one can choose a (5o > with 5o < 5 such that:there exists a smooth function 
g on with compact support such that g{Q, 0) = and D fl Bs^{xq) — {xo + aiBi + 
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(72 62 — sPxo I c""! + c^i + < ^0' > '^2)}, where Ci and are unit tangent vectors 
at Xq and orthogonal each other. 

One can choose a positive constant C in such a way that, for all a = {(Ti,a2) € 
it holds that \g{o')\ < C\a\. Let r > 1. We can easily sec that if s > C\a\, then 
s > \fTg{a I ^JT). This together with change of variables yields 

> [ e-'V\^'dads. 

Since 

lim / e-'V^^dads = [ e-'V^^dads < oo, 

T— ^oo Js>C\a\, |ct|2+s2<(V^5o)2 Js>C\a\ 

we have the desired conclusion. 

□ 



5.4 Claim 4 

Given 5 > one can choose a 5o > with Sq < 5 such that:there exists a smooth function 

h on with compact support such that h{0,0) — and dQ fl Bsq^ijo) = {ijq + (TiBi + 
(7262 — h{ai, o'2)i'yQ I (7^ + (t| + h{ai, 02)^ < Sq}, where 61 and 62 are unit tangent vectors 
at Ho and orthogonal each other. Note that h also satisfies that, for all ci = (cti, a2) G R^ 
\h{a)\ < C\a\, where C is a positive constant. We see that if |a| < 6o/{\/l + C^), then 
+ + /i((7i, (72)^ < Sq. This together with a change of variables yields 

/ e-^l2/o-2/l^c > f |V/i((7)Pcia 

J\a\<^5o 

Now one gets the desired conclusion. 

□ 
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